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Abstract 


We consider the problem of determining the free surface of a liquid in a capillary tube, 
and of a liquid drop, sitting first on a horizontal plane and then on more general surfaces. 
With some modifications, the method applies to the study of pendent drops and of rotating 
drops as well. 


Introduction 


Several capillarity phenomena, such as the rise of water in tubes of narrow bore, and the 
formation of liquid drops or bubbles, can conveniently be studied from the general point of 
view of the Calculus of Variations. Such a possibility, which clearly originates in the 
energy-minimizing character of the observed equilibrium configurations, has the remarkable 
advantage of providing a unified treatment of the mathematical questions that arise in a 
variety of particular phenomena. 

By using a well-known argument, based on the principle of virtual work, one is led to a 
variational formulation of the physical problem, in which a certain functional (representing 
the global energy of the system under consideration) has to be minimized, subject to some 
"natural" constraints, such as prescribed boundary conditions or fixed volume constraints. 

In general, the energy functional will consist of a "surface integral" plus a "volume inte- 
gral": the latter corresponds to body forces, of which gravity is a typical representative, 
while the former results, for example, from the consideration of the forces acting on the 
surface of separation between the liquid and the gas surrounding it. 

Now, the point is, that the classical definition of "surface area" is rather inadequate 
for treating this type of problem, mainly because it applies to smooth or Lipschitz-conti- 
nuous surfaces only - a class which is not closed under the usual limit operations. 

The difficulties arising from the presence of a surface integral become even more evi- 
dent when compared with the relatively simple treatment of the corresponding volume 
integral, which is generally well-defined on measurable sets and enjoys (at least in the 
simplest cases) nice variational properties. 

A satisfactory theory of surface area for a general class of surfaces of codimension one 
in ® n , n >. 2 , has been developed by E. De Giorgi in the fifties. 1 It is a remarkable fact 
that some classical questions, concerning the existence and regularity of capillary surfaces, 
have been answered only quite recently, using the variational techniques introduced by De 
Giorgi, or even more general methods pertaining to the field of Geometric Measure Theory. 2 

The definition of "functions of Bounded Variation", together with the main results of 
the corresponding BV-functions theory, will be recalled in the next section. As a first 
application of the theory, we shall discuss in section 2 the "standard” capillary problem, 
i.e. the determination of the free surface of a liquid in a thin tube of general cross 
section, which results from the simultaneous action of surface tension, boundary adhesion 
and gravity. It turns out that in this case the existence of the solution surface depends 
heavily on the validity of a simple geometric condition about the mean curvature of the 
boundary curve of the cross section of the capillary tube. Some particular examples of 
physical interest will also be discussed. 

Section 3 is devoted to the study of liquid drops, sitting on, or hanging from, a fixed 
horizontal plane. The symmetry of the solutions (which can actually be proved , as a conse- 
quence of a general symmetrization argument) now plays the chief role in deriving both the 
existence and the regularity of energy-minimizing configurations. When symmetry fails (this 
is the case, for example, when the "contact angle" between the drop and the plate is not 
constant, or when the supporting surface is not itself symmetric) , then more sophisticated 
methods must be used. Extensions in this direction will be outlined in section 4. 

We refer to the papers listed in the (fairly incomplete) bibliography at the end of the 
present paper for a deeper treatment of the subject, as well as for the discussion of 
related problems. 


Functions of bounded variation 

Given an open subset ft of IR n , n >. 2, we denote by BV(ft) the function space of Le- 
besgue integrable functions f over ft, whose distributional gradient Df is a vector 
measure with finite total variation on ft ; / ! Df | will denote the total variation of the 

measure Df , evaluated at A c ft . A 
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When A is open we obtain: 


1 1 Df | = sup(jf (x)divG(x)dx : G e Cq (A ; ]R n ) , | G | <.1 j 
A 


As a first result we can prove the lower semicontinuity of the map f ♦ / 1 Df j , with 
respect to the local convergence on ft ; that is, if fj «BV(ft) for every j and 

lim f | f . -f |dx = 0 

j-f + oo ' 

A 

for every A c; fi (i.e., A open and bounded with A c ft ) , then: 


j I Df | 

a 


£ lim inf j [ Df . | . 


Furthermore, we have the following compactness property: if fj eBV(ft) , and for every 
Ac eft and for every j it holds 

J | Df ^ | < c (A) 

A 

with c(A) independent of j , then there exists a subsequence of (f J , locally conver- 
ging in ft to some limit function f . 3 3 

If f f BV(P.) and ft has a Lipschitz-continuous boundary 3ft , then we can define the 
"trace" of f on 3ft (still denoted by f ), which is summable on 3ft and satisfies 4 

||f| < (1+L 2 ) 1/2 J | Df | + c (ft) 1 1 f i (1) 

3ft ft ft 


where L (the Lipschitz constant of 3ft ) and c(ft) depend only on the geometry of ft but 
not on f . 

By specializing the above definitions and properties to the case when f is the charac- 
teristic function $ E of a measurable set E c ft , we get a parallel theory of sets of 
finite perimeter, where by defin* *-ion: 


"perimeter of E 


in 


ft" 


= {|^ E I 

ft 


A straightforward application of the Gauss-Green theorem shows that this quantity coin- 
cides with the area of 3E " (1 , at least when 3E is a smooth (n-1 ) -dimensional surface 
in ft . The connection between BV-functions and sets of finite perimeter is given by the 
coarea formula: 


Df ! 



( 2 ) 


where f e BV(ft) and F t = (x t ft : f (x) < t) . 3 

Moreover 5 , if E = { (x, t) : x « ft, t < f (x) } , then itr holds 

J |d* e ! = j/i+[Df I 2 

ftx® ft 

where the second integral represents the total variation on ft of the vector measure, whose 
n + 1 components are respectively the Lebesgue measure on K n , and the distribution 
derivatives D^f, i=1,...,n . When f is Lipschitz-continuous on ft , this yields of 
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course the area of the graph of f over S 1 . 

Capillary surfaces In cylindrical vertical tubes 

Let us consider a capillary tube, open at both ends and partially immersed into a liquid; 
for simplicity, we may assume that the liquid rises in the tube, so that - in view of the 
preceding discussion - the energy corresponding to a certain configuration of the liquid 
within the tube, described by the graph a function f > 0 , can be expressed as 


*J(f> 



(3) 


Here, £2 (open and bounded in IR 2 , with Llpschitz boundary 3ft ) denotes the cross 

section of the tube, f t iiV(il) , and k,v are physical constants, with < N 0 and 0 < v < 1 . 

It is easily seen that f 0 is the trivial solution of *2(*) -* min , when v < 0 , while 

inf ‘'U • } = when v ' 1 . 

We can immediately check that a configuration of minimal energy (i.e., a solution to the 
problem *"{(■) -* min ) satisfies the equations 


( div Tf(x,y) = kf(x,y) in si 
l Tf(x,y) • n(x,y) = v on 3d 


where 


Tf = 


Df 


1 + 


and n(x,y) denotes the outward unit normal to 3s2 at (x,y) 


,1 


C 1 (IS) 


This is true if, for example, 'C is of class C and f t C (Si) n 

This way, we realize that the moan curvature of a capillary surface is, at any point 
(x,y,i (x,y) ) t si'®, proportional (with constant k ) to its height above the reference 
plane z = 0 , and that v corresponds to the cosine of the contact angle between the sur- 
face and the walls of the tube. 


The classical approach consisted in solving the above system of equations in the special 
case when s'; was a disc of radius R , and f = f(r), r = (x 2 + y 2 ) '/ 2 (axially symmetric 
solutions) . In this case, one is led to the ordinary differential equation 



k r f ( r ) 


for 0 -- r s R 


with the boundary conditions: f r (0) = 0, f r (R) = v- (1-v 2 ) ~1/2 . 

In order to prove the existence of a solution f to the problem 'J ( • ) -* min , we make 
the following assumption on the domain si : f 



3 si 





(4) 


for every Eos;, with a i 1 and 6 s 0 . From (4) and the coarea formula (2) we get 

jf v aj i Df I+ d jf 


for every f 6 BV(sl) with f 1 0 . By introducing this last inequality in (3) we find: 

2 2 

-;f(f5 > (1-va) i Df | + ^jf 2 - — - - In’; , 

Si si 


where Id! denotes the Lebesgue measure of 
inf }(■) 2 - v 2 6 2 I si ! / k , while if v » 1/a , 
ing ' J(fj) -* inf r J( ■ ) , we obtain: 


si . In particular, if v 1/a 
then for any minimizing sequence 


we obtain 
fj , satisfy- 


|Df j! 


const . 
1 - vi 



»: 


const . 
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so that, in view of the results in section 1, we conclude that a subsequence of { f -j } con- 
verges to a function f , which clearly solves our problem; as a consequence of the strict 

convexity of the energy functional, such a solution is in fact unique. 

The hypothesis (4) now comes into discussion: clearly, it is implied by the trace esti- 
mate (1), so that capillary surfaces always exist for v in the range 0 < v < (1+L : ) -1 ^ , 
with L = Lipschitz constant of 90 ; k that this is in fact* an "almost necessary" condition 
can easily be seen with the aid of some simple examples. 7 For a domain 0 in the form of 
a circular sector, no solution with bounded energy can exist when v > ( 1 +L ~ ^ (that is, 
when 0 + 2y < tt , where P is the angle of the sector and > is the contact angle) . The 
discontinuity at 0 + 2> = it is also confirmed by physical experiments. 7 

On the other hand, 6 interior corners do not affect the solution: in the special case when 
d satisfies an Internal Sphere Condition of radius R N 0 (i.e., when each point in 0 
belongs to some ball Br c si ) , then (4) holds with oc = 1 , and it can be shown, by using 
a result of Concus and Finn, 6 that the solution exists for every v with 0 < « < 1 . 

In conclusion, we remark that the variational method does not work when x =* 0 , i.e. in 
the absence of gravity; one can actually show* that the solutions f k , corresponding to 
values x 0 , go uniformly to +°° in i? , as x •» 0 + . 

Sessile drops and pendent drops 

The energy of a liquid drop, sitting on the horizontal plane (z=0) in 1R 5 , can be 
written in the following way: 


(E) 



z^O 


+ v 


*E * 


z=0 



z s 0 


(5) 


where E denotes the region of the half-space <z > 0) occupied by the liquid. The first 
integral in (5) represents the area of the free boundary of the drop, the second integral 
gives the area of the region of contact, and the third integral corresponds to gravity. As 
usual, we assume x . N 0 and v t ( — 1,13 , since for v = -1 no solution can occur. The 
same functicnal, with x 0 , represents the energy of a pendent drop; in both cases, a 
volume constraint has to be imposed, namely |e| = v ' 0 . 

Now, it can be shown 18 that by replacing each horizontal section of a given configuration 
E by a disc of the same area, centered on the z-axis, a new configuration E 1 results, 
which is of less energy than E (in fact, \ (E) = ^.(E 1 ) iff E is already axially 

symmetric) . 

From this fact, by using the obvious estimate 



z=0 z^O 


which implies 


(E) 1 


1+v 




for every E , if k 1 0 , one gets easily the existence of a minimum of 'Jx ( ‘ ) , when 

x > 0 . On the other hand, when x < o the situation is completely different, and we can 
look only for local solutions of '7x ( • ) ‘ min , since clearly inf '~1 K ( • ) = in this case. 

To this aim, we introduce the following definition: E is a local minimum of the energy 

functional (a pendent drop) if ! E I = v , and there exist T N 0 an3 u f ( 0 . 1 ) such that' 

E is contained in the region to < z < «t) and, for every F c (0 < z < T> with I F I « v , 

there holds ‘^ k (E) £ "? k -(F) • 

In order to prove the existence of pendent drops, we start by observing that, when x = 0 
(that is, in the absence of gravity), the minimum E c of ‘ ? 0 ( ’ ) is a portion of a ball, 
completely determined from the data v and v . With this in mind, we choose T 
than the maximum height of E 0 and find the solution 
restricted to the configurations E c (o < z <• T^ . 

An easy calculation shows that Ex ■* E 0 as x -* o 
results concerning the convergence of surfaces of prescribed mean curvature, we get that for 
I x | small enough, the solution Ex is actually contained in (0 v z x ot) for a suitable 
a t(0,1), thus concluding the proof of the existence of pendent drop^ of given volume, in 
a weak gravitational field. Estimates on the effective smallness of Im 
tly known. 

For example, Giusti 1- showed that pendent drops exist, if the product |x|v r/ ’ does not 
exceed a constant, which can be written down expl icitly, and which depends only on the value 
of the contact angle between the free surface of the drop and the plane { z=0 ) . 


(•) 


greater 
min , 


of the problem 
as a consequence of general 


are also explici- 
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Concluding comments 


The regularity (analyticity) of the equilibrium surface of a liquid in a capillary tube 
can be derived from general regularity results for hyper surfaces of least area or, more ge- 
nerally, of prescribed mean curvature in ]R n . 13 The boundary behaviour of the solutions 
has also been studied, and several results in this direction are presently known. 14 

As far as the regularity of liquid drops is concerned, we note firstly that the free 
surface of a drop sitting on a horizontal plane, being rotationally symmetric, can also be 
described (locally) as the graph of a suitable function f , defined on a 2-dimensional 
domain. It turns out that f minimizes a functional of the type: 


gravity + X 


where the Lagrange multiplier * takes into account the volume constraint. From this, the 
interior analyticity of the free surface of the drop follows at once. 

Secondly, it is not difficult to prove that the configuration representing a sessile drop 
is a convex set. 6 We then conclude that its free surface is in fact smooth up to the plane 
{ z = 0 } , and that the cosine of the contact angle coincides with . 

Much more difficult is the study of liquid drops, when symmetry ceases to exist. However 
some partial results have recently been obtained, such as the existence and regularity of 
solutions, corresponding to drops sitting q>n Jor hanging from) a surface z = >i(x,y) in 
1R 3 , satisfying C (r) +“’ for r = (x'+y 1 ) 1 -*■ +“ . From this "growth condition at infini- 

ty", which however is not satisfied in a number of interesting and still open situations, 
the necessary compactness results can easily be derived. 17 

A theorem, which extends to minima of area-like functionals, subject to a volume con- 
straint, the regularity theory for minimal boundaries, is now available. 17 ' 18 A technique, 
developed in connection with its proof, can also be used to prove the existence of liquid 
drops, in rotation around an axis through the center of mass, and held together by surface 
tension. 

In its simplest formulation, the problem asks for a local minimum of the functional 


■ 10 


(E) 



ul 

) 

1R 


(x 


+y ‘ ) dxdydz 


U) > 0 


( 6 ) 


subject to a volume constraint (|e| = 1) and to a further constraint about the center of 
mass (which must coincide with the origin of the space) . 

The presence of the kinetic energy excludes, in general, the symmetry of ■ :ie solutions. 
Anyway, following the treatment of the pendent drop problem, one can show 18 the existence of 
relative minima of the energy functional, when u is small enough, that is, when the rota- 
tion is sufficiently slow. 

The proof of this result, which can only be outlined here, proceeds as follows: firstly, 
we define E to be a local minimum of the energy functional (6) , i f there exists 
R ■> R 0 = ( 3/4 tt ) 1/3 such that E cc c R and, for every admissible F cc Cr , satisfying the 
above constraints, it holds i_ ^(F) ■ 

Here, Cr denotes the cylindrical container 


C R = | (x , y , z) tr 3R 3 : x ; +y : £R', |z| £ R) 


Then, we can easily prove that for every u > 0 , and every R > R 0 , there exists a solu- 
tion E^ to the problem '? U (E) m in > restricted to the admissible configurations E c C R 
Such a solution, however, will not generally satisfy the condition E^ cc c R for a local 
minimum. 

Next, we observe that for a fixed R > R 0 , the solutions E u we have found in this way, 
converge, as o -*• 0 + , to the ball E 0 = {(x,y,z) : x‘+y ; +z : < Rq) cc C r . The local conver- 
gence E u - E 0 is unfortunately too weak to conclude that E w cc c R for oi small enough. 
But we can prove, and this is actually the crucial point of the entire demonstration, that 
if uj is sufficiently small (and positive), then there exists a value r <-(R 0 ,R) such that 


?c r 

From this result, the fact that E^ccCr for w small can be proved as follows. Define: 
G = E^ n Cj- and a = [E^ - C r | 
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and choose F = (1-a)' 1/3 G (that is, F is constructed by expanding G , with a coefficient 
(1-a) -1 ^ > 1 ). Clearly, | F j = 1 . Since a •+ 0 as w •* 0 , it is clear that when w is 
small, by means of a suitable translation F ■+ F we obtain eventually an admissible confi- 
guration ? c Cr , whence 

7 u (Ej<7 u (f) (7) 

On the other hand, the difference ^(E^) - ( F ) can easily be estimated from below by the 

quantity: 

[l - (1 - a)' 2/3 ] j|D<t> G l + const. a + const. 1 a 2/3 
TR 3 

Now, from the Taylor expansion of the coefficient between square brackets, we conclude that 
r % i (E u ) - J^{F) >0 if a is positive and small, which would contradict (7). Therefore, 
for ui small enough, a must be 0 , so that the corresponding solution is contained 

in Cj- Cp . 

In conclusion, we remark that the preceding argument, when used in connection with gene- 
ral regularization techniques, 20 allows the proof of the analyticity of the solution E u as 
well . 
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